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Plasma  Irregularities  in  the  D-Region  Ionosphere  in  Asso¬ 
ciation  with  Sprite  Streamer  Initiation 

Jhmqj  Qin'%  Victor  P,  Pasko1*  Matthew  G.  McHarg2  and  Hans  C.  Stenbaek-Nielsen3 

1  Com  munications  and  Space  Sciences  Laboratory,  Department  of  Electrical  Engineering,  Penn 
State  University ,  University  Park,  Pennsyl vania,  USA 

2 Department  of  Physics,  United  States  Air  Force  Academy,  Co/orado  Spring  Colorado,  USA 
Geophysical  Institute,  University  of  Alaska  Fairbanks,  Fairbanks,  Alaska,  USA 

Sprites  are  spectacular  optical  emissions  in  the  mesosphere  i&ctuced  by  transient  lightning 
electric  Adds  above  thunderstorms.  Although  the  streamer  nature  of  sprites  lias  been  gener¬ 
ally  accepted,  how  these  filamentary  plasmas  are  initiated  in  the  lower  ionosphere  remains  a 
subject  of  active  research  and  existing  theories  differ  fundamentally  from  each  other.  Here 
wc  present  a  combination  of  observational  and  modeling  results  showing  solid  evidence  of 
pre-existing  plasma  irregularities  in  association  with  sprite  streamer  initiation  In  the  l>- 
region  ionosphere.  The  high-speed  video  observations  show  that  prior  to  streamer  initia¬ 
tion,  spatial  structures  discernible  in  overall  diffuse  emissions  of  sprite  halo  descend  rapidly 
with  the  halo,  but  slow  down  and  stop  to  form  the  stationary  glow  in  the  vicinity  of  the 
streamer  onset,  from  where  streamers  suddenly  emerge.  The  modeling  results  reproduce  the 
above-mentioned  halo  dynamics  and  demonstrate  that  the  descending  structures  in  the  sprite 
halos  are  the  optical  manifestations  of  preexisting  plasma  irregularities  in  the  D- region  iono¬ 
sphere,  which  are  necessary  for  streamer  inception.  An  image  processing  algorithm  is  devel¬ 


oped  to  reconstruct  the  plasma  irregularities,  showing  that  these  plasma  irregularities  have  a 
vertical  dimension  of  a  few  kilometers  and  a  horizontal  dimension  of  a  few  hundred  meters. 
The  plasma  irregularities  might  have  been  produced  by  thunderstorm  or  meteor  effects  on 
the  l)-region  ionosphere. 

From  the  Earth's  surface  to  ils  ionosphere,  the  planet  s  atmosphere  experiences  an  expc^ 
nemial  decrease  of  its  air  density  and  an  exponential  increase  of  its  electron  density,  leading  to  a 
transition  region  with  unique  properties  in  die  altitude  range  of  60  io  90  km,  known  as  the  D  region 
of  die  ionosphere,  between  the  planet’s  lower  atmosphere  with  dense  ait  and  its  highly  conducting 
ionosphere.  For  tropospheric  thunderstorms,  low  electron  density  in  this  transition  region  enables 
the  penetration  of  lightning  electric  fields,  lead i tig  to  many  not  yet  well  understood  thunderstorm 
effects  such  as  electron  density  fluctuations  and  depletions3’2,  as  well  as  occasionally  the  formal  ion 
of  spectacular  sprite  discharges1.  Moreover,  for  meteor  events,  the  D  region  is  usually  Ihe  lowest 
altitude  region  that  meteoroids  can  reach  due  to  significant  friction  caused  by  the  rapid  increase  of 
air  density,  resulting  in  plasma  irregularities  observed  as  low-altitude  meteor  nail  echoes  in  radar 
observations4. 

Predicted  in  1925  and  first  documented  in  19E9,  sprite  discharges  originating  from  the  D 
region  ionosphere  have  been  an  active  research  area  in  the  pasi  two  decades5-41.  Early  sprite  theory 
emphasized  the  importance  of  thunderstorm  activity  since  it  was  realized  soon  after  the  discovery 
I  hat  intense  cloud-to-ground  lightning  discharges  associated  with  large  charge  moment  changes 
(j,e„  chaige  times  altitude  from  which  it  was  removed)  are  necessary  for  their  production7.  Fol- 
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lowing  that,  the  initiation  of  sprite  streamers  was  interpreted  in  a  modeling  study  as  a  result  of  the 
sharpening  and  collapse  of  the  screening- ionization  wave  associated  with  the  sprite  halo  when  the 
lightning  charge  moment  change  was  as  large  as  1300  C  km15.  That  theory  encounters  difficulties 
in  explaining  the  temporal  and  spatial  relations  between  sprite  halo  and  sprite  streamers^*  and  ex¬ 
plaining  them  using  experimentally  observed  (relatively  small)  charge  moment  changes  with  an 
empirical  minimum  of  300  C  km|r>.  A  different  theory  proposed  recently  overcomes  these  difficul¬ 
ties  and  indicates  that  the  presence  of  strong  plasma  iiTegularities/inhomogeneities  in  the  D  region 
ionosphere  is  a  necessary  condition  for  the  initiation  of  sprite  streamers5"*  in  agreement  with 
recent  experimental  findings  showing  that  a  large  charge  moment  change  is  a  necessary  but  not 
-sufficient  condition  for  sprites1’.  However,  up  to  dale*  no  solid  evidence  of  D-region  plasma  irreg¬ 
ularities  in  association  w  ith  sprite  streamer  initiation  has  been  presented*  and  if  these  irregularities 
exist  Iheir  characteristics  and  origin(s)  remain  unknown*  In  this  study,  we  analyze  high-speed  video 
records  of  sprite-halo  events  and  focus  specifically  on  the  halo  dynamics  prior  to  sprite  streamer 
initiation-  By  comparing  observational  and  numerical  results,  we  demonstrate  that  the  descend¬ 
ing  halo  structures  recorded  in  the  high-speed  video  observations,  from  which  sprite  streamers  are 
initiated,  are  optical  manifestation  of  Site  pre-existing  plasma  irregularities.  An  image  processing 
algorithm  is  then  developed  to  reconstruct  the  characteristic  shape  and  size  of  those  plasma  irreg¬ 
ularities*  showing  that  sprite  halos  can  be  effectively  used  to  probe  the  plasma  irregularities  in  the 
lower  ionosphere.  We  further  discuss  origins  of  these  plasma  irregularities,  suggesting  possible 
relations  between  sprites  and  other  thunderstorm  effects  or  meteor  events. 


1  Origin  of  (he  Descending  Halo  Structures 

In  video  observations*  sprite-halo  events  exhibit  a  brief  descending  diffuse  glow  in  the  shape  of 
a  pancake  with  diameters  up  to  ^80  km  near  ^7 5  km  altitude*  referred  to  as  a  sprite  halo13,  and 
develop  into  fine-structure  filaments  with  diameters  up  to  several  hundred  meters  in  the  altitude 
range  of —40  to  ^90  km,  commonly  referred  to  as  sprite  streamers1'3.  Recent  imaging  of  sprites  at 
16,000  fp*  indicates  that  prior  to  streamer  initiation*  spatial  structures  in  halo  emissions  descend 
rapidly  with  the  sprite  halo*  Ixu  slow  down  and  stop  to  form  the  stationary  glowr  in  the  vicinity 
of  the  streamer  onset,  from  where  streamers  suddenly  emerge  and  then  accelerate  exponentially3. 
The  above-mentioned  halo  dynamics  is  illustrated  in  Figure  1. 

The  descending  halo  structures  appear  to  be  of  esse  alia  I  importance  for  the  initiation  of  sprite 
streamers.  To  understand  the  origin  of  those  structures*  it  should  be  first  emphasized  that  in  the 
transient  process  of  halo  development  on  a  typical  time  scale  of  M  ms.  free  electrons  in  the  lower 
ionosphere  move  only  —100  m  assuming  an  electron  drift  velocity  of  I0r>  m/s,  which  is  negligible 
when  compared  to  the  descending  distance  of  halo  structures  on  the  order  of  -10  km.  In  particular, 
in  all  halo  events  presented  in  this  work,  the  electrons  move  upwards  (i.e,  Jri  an  opposite  direction 
to  the  downward  movement  of  llic  halos)  due  to  the  positive  polarity  of  the  causative  lightning 
dischaiges.  The  observed  descent  of  a  sprite  halo  as  an  entire  object  is  due  to  the  increase  of  the 
large-scale  lightning  electric  field  at  low  altitudes  and  fast  relaxation  at  high  altitudes  during  (he 
causative  lightning  discharge,  leading  to  the  lowering  of  the  thin  altitude  region  of  large  reduced 
electric  field  EfN,  where  N  is  the  air  density  at  the  attitude  of  interest  (he,*  the  large-scale  halo  is 
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optical  manifestation  of  the  region  of  large  reduced  electric  field  in  the  lower  ionosphere}9.  Note 
that  large  reduced  electric  held  leads  to  high  rate  of  electron  impact  excitation  of  Nir  which  is  the 
main  physical  process  accounting  for  the  halo  optical  emissions.  Thus  the  observed  descent  of 
halo  structures,  proceeding  at  the  same  pace  as  the  large-scale  halo  emissions,  is  also  due  to  the 
lowering  of  the  thin  altitude  region  oflarge  reduced  electric  field.  I  ti  other  words,  at  each  moment 
of  time,  the  vertical  dimension  of  the  halo  structures  should  be  defined  by  the  vertical  structuring 
of  the  lightning  reduced  electric  field.  However,  unlike  its  rapid  spatial  variation  in  the  vertical 
direction,  [he  lightning  electric  field  varies  slowly  on  a  spatial  scale  of  tens  of  kilometers  in  the 
horizontal  direction,  which  cannot  lead  to  the  observed  small-scale  horizontal  fluctuations  in  the 
halo  emissions.  Since  the  electron  impact  excitation  rate  is  determined  by  the  reduced  electric 
field  along  with  the  local  electron  density,  we  propose  that  the  physical  parameter  that  determines 
the  horizontal  dimension  of  the  halo  structures  should  be  the  fluctuations  in  the  ambient  electron 
density  (i.e.?  ambient  plasma  irregularities}. 

To  demonstrate  Ehat  the  downward  progression  of  the  high  reduced  field  region  and  the  pres¬ 
ence  of  plasma  irregularities  can  account  for  the  observed  descending  halo  structures,  we  use  a 
two-dimensional  cylindrical!}'  symmetric  plasma  fluid  model  (see  Methods  section)  to  numerically 
reproduce  the  dynamics  of  the  descending  halo  structures*  Figure  2a  shows  the  ambient  electron 
density  and  the  plasma  irregularity  assumed  in  the  simulation  domain,  lire  plasma  irregularity  has 
a  Gaussian  distribution  in  the  horizontal  direction  with  a  peak  value  of  IQQcm^  aEid  a  characteris¬ 
tic  size  of  2  km.  In  the  vertical  direction,  Ehe  plasma  irregularity  is  homogeneous  and  has  a  length 
of —  12  km,  comparable  to  the  deseed  mg  distance  of  halo  structures  shown  in  Figure  I.  Panels  b 


and  c  in  Figure  2  show  the  sprite  halo  emissions  from  the  hirst  positive  band  system  of  at  £=0.7 
ms  and  1.2  ms,  respectively.  It  is  dear  that  at  each  moment  of  rime,  only  a  small  portion  of  the 
plasma  irregularity  produces  emissions  that  are  brighter  than  those  originating  from  other  regions 
of  the  halo.  The  horizontal  dimension  of  this  portion  is  approximately  equal  to  the  assumed  hori¬ 
zontal  size  of  the  plasma  irregularity.  Examination  of  the  spatial  variaiion  of  the  lightning  electric 
field  at  each  moment  of  time  shows  that  the  bright  portion  of  the  plasma  irregularity  is  located  at 
the  altitude  region  with  the  largest  lightning  reduced  electric  field,  and  the  vertical  dimension  of  the 
bright  portion  is  determined  by  ihe  vertical  size  of  this  altitude  region.  Figure  2d  shows  a  1,4-ms 
image  time  series  of  the  modeled  halo  structure.  We  can  see  dearly  that  the  bright  potion  of  the 
plasma  irregularity  descends  rapidly  at  tire  same  pace  as  the  remaining  parts  of  halo  emissions,  but 
the  descent  slows  down  due  to  the  decrease  of  lightning  current,  and  at  f=l*4  ms  sprite  streamers 
stall  forming.  The  speed  of  the  desccni  (7,6x  10°  m/s>  is  comparable  to  that  of  the  halo  structure 
(6,4 x  1 015  m/s)  shown  in  Figure  1,  We  note  that  tire  streamer  starts  after  the  halo  has  faded  due  to 
a  longer  time  scale  (>  1  ms)  required  for  streamer  initiation  at  -75  km  altitude  when  compared  to 
the  duration  of  the  halo  (-1  ms  for  typical  lightning  return  strokes.)  It  also  needs  to  be  emphasized 
that  although  the  model  indicates  the  inception  of  a  sprite  streamer  at  id  1.4  ms,  it  cannot  accu¬ 
rately  shnulate  the  dynamics  of  the  sprite  streamer  at  later  moments  of  time  due  to  the  low  spatial 
resolution  1 20  m)  used  in  the  large  simulation  domain*  In  the  head  of  a  streamer  at  —70  km 
altitude,  the  electros*  density  could  change  one  order  of  magnitude  on  a  spatial  scale  of  —2  m.  This 
high  electron  density  gradient  requires  sub-meter  spatial  resolution  to  simulate,  which  is  beyond  a 
practical  resolution  possible  in  a  large  simulation  domain  used.  Nevertheless,  as  shown  jn  Figure  2, 
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with  a  spatial  resolution  of  M20  m,  the  model  is  fully  capable  of  reproducing  large-scale  plasma 
dynamics  that  involves  only  low  electron  density  gradients,  such  as  the  experimentally  observed 
descending  halo  structures.  Based  on  Else  above  analysis  and  comparison*  it  can  be  concluded  that 
the  halo  structures  are  originated  from  preexisting  D-region  plasma  irregularities*  and  the  descent 
of  the  halo  structures  is  due  to  the  lowering  of  the  thin  altitude  region  oflarge  reduced  electric  field 
during  (he  charge  removal  process  of  the  causative  lightning  discharge. 

2  Reconstruction  of  the  Pre-existing  Plasma  Irregularities 

Since  at  different  moments  of  time  during  (he  lightning  discharge,  different  portions  of  the  plasma 
Irregularity  are  effectively  “illuminated"  by  the  lightning  quasi-, static  electric  field,  the  descending 
halo  structures  recorded  by  high-speed  video  cameras  can  be  used  to  reconstruct  the  shape  and  size 
of  (he  plasma  irregularities  in  the  lower  ionosphere.  Figure  3a  shows  the  plasma  irregularity  in  the 
modeled  sprite-halo  event  reconstructed  using  image  series  from  t=0  ms  to  1*4  ms  (see  Figure 
2d)  and  the  image-processing  algorithm  described  in  the  Methods  section.  Comparison  between 
Figures  2a  anti  3a  indicates  that  the  plasma  irregularity  in  (he  simulation  domain  can  be  well 
reconstructed.  Figure  3b  shows  the  D-region  plasma  irregularity  reconstructed  using  high-speed 
video  records  of  the  sprite-halo  event  shown  in  Figure  I.  We  emphasize  that  the  reconstructed 
image  does  not  contain  information  about  the  vertical  density  variation,  but  only  shows  the  shape 
and  size  of  the  plasma  irregularity  as  manifested  by  optical  emissions.  Nevertheless*  important 
information  can  be  derived  from  it  First,  sprite  streamers  can  be  easily  Identified  because  they 
are  much  brighter  than  the  plasma  irregularity  and  the  halo.  Second,  the  right  streamer  is  initiated 
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from  the  lower  tip  of  the  plasma  irregularity.  Plasma  irregularity  responsible  for  the  left  streamer 
cannot  be  reconstructed  because  if  presents  it  should  be  located  above  the  available  field  of  view 
of  the  camera.  Third,  the  size  of  the  plasma  irregularity  can  be  estimated  using  two  different 
approaches,  ( i)  Based  on  (he  altitude  range  of  the  entire  image  which  is  calculated  assuming  the 
sprite  streamer  was  directly  above  the  causative  lightning  discharge,  the  vertical  and  horizontal 
size  of  the  plasma  irregularity  are,  respectively,  ^5  km  and  ^  1  km.  (2)  Based  on  a  100  m  typical 
diameter  of  sprite  streamers  documented  in  the  existing  literal ure15,  the  vertical  and  horizontal  size 
of  the  plasma  irregularity  are*  respectively,  ~2.5  km  and  -^0.5  km.  These  two  estimates  are  roughly 
consistent  with  each  other.  Note  that  it  is  possible  that  the  plasma  irregularity  has  a  significantly 
longer  vertical  dimension  than  that  revealed  in  Figure  3b*  because  the  lightning  electric  field  cannot 
penetrate  into  higher  altitudes  so  that  even  if  a  significant  portion  exists  at  high  altitudes*  it  cannot 
be  “illuminated”* 

To  further  demonstrate  the  characteristics  of  plasma  irregularities  in  association  with  sprite 
streamer  initiation*  we  analyze  another  sprite-halo  event  observed  on  20  July  2012*  in  which  multi¬ 
ple  plasma  irregularities  can  be  identified  in  the  halo.  For  this  event*  the  halo  emissions  first  appear 
at  06:27:17:157  860  UT  (i=0  ms).  Figure  4a  shows  the  halo  emissions  at  f=L04  ms.  Five  bright 
halo  structures  are  identifiable  at  this  moment  of  time,  and  in  the  video  record  they  all  experience 
descending  motion  similar  to  that  shown  in  Figure  Id.  A I  £=1,28  ms,  a  sprite  streamer  was  already 
initiated  from  the  leftmost  halo  strucfure.  For  all  (he  olher  four  halo  structures,  £=1*28  ms  is  the 
moment  just  before  the  initiation  of  sprite  streamers.  Figure  4b  shows  (lie  five  plasma  irregular¬ 
ities  reconstructed  using  all  video  frames  from  t  -  0  to  1.28  ms.  These  plasma  irregularities  are 
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respectively  responsible  for  the  live  descending  halo  structures  (see  Figure  4a)  and  the  follow¬ 
ing  live  sprites  (not  all  shown  In  Figure  4b),  As  an  additional  evidence.  Figure  4c  shows  a  halo 
event  w  ith  no  following  sprite  streamers.  This  halo  event  was  observed  on  14  July  2010  starting  at 
08:32:3 1.594  409  LIT,  and  it  has  a  typical  pancake  shape,  which  is  defined  by  the  spatial  variation 
of  the  large-scale  lightning  electric  field.  In  the  reconstructed  image  shown  in  Figure  4d,  no  plasma 
irregularities  and  sprite  streamers  can  be  identified*  which  is  an  indirect  evidence  suggesting  the 
importance  of  plasma  irregularities  for  the  initiation  of  sprite  streamers. 

A  recently  proposed  theory  suggests  that  the  presence  of  strong  plasma  irregularities  in  the 
D-region  ionosphere  is  a  necessary  condition  for  the  initiation  of  sprite  streamers'^ 1 1  f  and  modeling 
results  show  that  vertically  elOFigaled  plasma  columns  are  more  favorable  for  streamer  initiation 
when  compared  to  spherical  irregularities1^17.  The  above  reconstructed  images  demonst rate  for  the 
lirst  time  solid  evidence  of  kilometer-scale  plasma  irregularities  in  association  with  sprite  streamer 
initiation  in  the  D- region  ionosphere*  However,  it  is  important  to  emphasize  that  not  all  the  plasma 
irregularities  in  association  with  sprite  streamer  initiation  can  be  observed  as  descending  halo 
structures  in  video  records  of  sprite  events.  Two  reasons  can  lead  to  the  absence  of  descending 
halo  structures  in  sprite  observations,  Firsi,  the  descending  halo  structures  emphasized  in  the 
present  work  are  due  to  a  significant  vertical  dimension  of  the  plasma  irregularities  leading  (o  the 
effective  illumination  of  a  small  portion  of  the  plasma  irregularities  at  each  moment  of  time  by  die 
lightning  electric  field.  In  the  case  when  the  plasma  irregularities  are  approximately  spherical  and 
have  a  compact  size  of  tens  to  hundreds  of  meters,  which  are  capable  of  initiating  sprite  streamers 
according  to  Lite  previous  modeling  studies1118,  the  plasma  irregularities  could  only  be  recorded 


as  a  local  enhancement  of  the  halo  luminosity  around  the  origin  of  sprite  streamers.  Second,  the 
brightness  of  halos  depends  on  the  impulsiveness  of  the  causative  lightning  return  strokes1*.  In 
some  sprite  events,  [he  halo  and  halo  structures  can  remain  sub- visual  if  the  return  stroke  is  not 
impulsive  enough.  We  note  that  a  relatively  bright  but  not  saturated  large-scale  halo  is  ideal  for  the 
reconstruction  of  the  plasma  irregularities*  H  has  also  been  verified  that  in  some  oases  when  only 
the  descending  halo  structures  are  visible,  the  reconstruction  is  also  possible  since  halo  structures 
are  brighter  lhan  the  other  parts  of  die  halo,  that  remain  sub-visual. 

3  Possi  blc  Or  i  gins  of  the  Reconst  meted  Plasma  Irregulari  ti  cs 

According  to  previous  modeling  studies1618*  the  electron  density  in  plasma  irregularities  needs 
to  be  ^2  to  5  orders  of  magnitude  higher  than  the  ambient  electron  density  in  order  to  initiate 
sprite  streamers,  depending  on  the  onset  altitude  of  the  sprite  streamers  (at  lower  altitudes,  the 
difference  needs  to  be  larger),  St  is  reasonable  to  believe  that  such  strong  perturbations  can  not 
persist  in  a  quiet  lower  ionosphere*  and  it  is  important  to  emphasize  that  the  origin(s)  of  Ihose 
plasma  irregularities  cannot  be  determined  using  the  results  presented  in  this  work.  The  several 
possible  origins  are  discussed  below, 

(1)  Thunderstorm  effects.  Tropospheric  thunderstorms  are  well-known  to  be  able  to  disturb 
the  lower  ionosphere  by  convective  atmospheric  gravity  waves  and  by  lightning  electric  fields, 
leading  to  conductivity  enhancement,  electron  density  depletions*  and  D-layer  re  fleet  ion  height 
splitting1,2’7'19.  However,  most  of  those  are  large  spatial  scale  processes.  The  only  known  ill  unde  r- 
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storm  effect  that  could  produce  the  kilometer-scale  plasma  irregularities  with  a  density  enhance¬ 
ment  of  at  least  —2  orders  of  magnitude  is  sprites  themselves.  Previous  modeling  results  suggest 
that  the  electron  density  enhancement  in  sprite  bodies  cat)  typically  last  for  approximately  hundred 
seconds  at  —75  km  altitude™.  The  remnant  of  a  previous  sprite  body  could  indeed  facilitate  the  ini¬ 
tiation  of  sprites  at  (eus  of  milliseconds  later  according  to  video  observations21,22.  To  examine  the 
possibility  of  previous  sprites  producing  the  plasma  irregularities  observed  in  Figures  J  and  4a,  we 
have  analyzed  the  images  of  sprites  captured  prior  to  those  events  on  the  same  nights.  We  find  that 
no  previous  sprites  occurred  at  the  same  location  as  the  one  shown  in  Figure  l*  For  the  event  shown 
in  Figure  4a,  a  sprite  event  occurred  —  ID  minutes  before,  and  the  bright  sprite  columns  captured  by 
the  same  camera  largely  overlap  with  the  sprites  shown  in  Figure  4a,  Unfortunately,  images  cap¬ 
tured  using  another  camera  missed  most  parts  of  the  two  events  so  that  the  exact  coincidence  of  the 
sprite  locations  cannot  be  verified  using  a  tri angulation  technique.  The  -10  mi smies  time  separa¬ 
tion  is  much  longer  than  the  —  100  s  relaxation  time  of  the  electron  density  enhancement  at  —75  km 
altitude  estimated  in  modeling  of  sprite  chemistry20.  On  the  other  hand,  recent  coordinated  radio 
and  optical  observations  show  that  long  recovery  early  VLF  events,  caused  by  long* lasting  (a  few 
hundred  seconds  to  tens  of  minutes)  D-region  ionospheric  conductivity  modifications,  occurred 
often  coi nci dent ly  with  bright  sprites23,  It  has  been  suggested  (hat  the  long- lasting  conductivity 
modifications  may  have  been  produced  by  those  concurrent  sprites23.  Although  how  the  electron 
density  enhancement  in  sprites  could  last  so  long  is  not  dear,  the  plasma  irregularities  shown  in 
Figure  4a  might  be  related  to  the  previous  sprite  event  that  occurred  -10  minutes  before. 


(2)  Meteor  events.  Considering  the  requirement  of  several  orders  of  magnitude  plasma  den¬ 


sity  enhancement  for  streamer  initiation,  those  plasma  irregularities  might  have  been  produced 
by  meteor  events.  Meteors  had  been  suggested  to  be  a  possible  factor  for  the  initiation  of  sprite 
streamers  in  early  sprite  resea ncht  but  soon  the  idea  was  discarded  because  using  a  visually  ob¬ 
served  meteor  flux  density  of  8x  10“ 7  km_as_1  it  wax  estimated  that  for  a  sprite  with  a  horizontal 
cross-sectional  area  of  —175  km2  one  would  expect  only  —ID"1  sporadic  meteors  io  intersect  a 
sprite  volume  in  a  1  s  interval  prior  to  the  optical  observation  of  the  sprite24.  However,  radar  obser¬ 
vations  show  that  visual  and  sub- visual  meteors  are  continuously  entering  the  Barth's  upper  atmo¬ 
sphere  with  a  much  larger  (lux  density  estimated  to  be  -IG“*  km~  V]  globally2*.  Moreover,  the 
entire  volume  of  a  large-scale  halo,  which  covers  a  horizontal  cross-sectional  area  of  5000 

km2  assuming  a  radius  7?ilrtio  of  40  km*  should  be  used  in  the  estimate  because  sprites  could  have 
lateral  offsets  as  large  as  —40-50  km  with  inspect  to  the  causative  lightning  discharge26.  With 
these  different  values,  it  is  estimated  that  on  average  0.05  meteors  could  Intersect  a  large-scale 
halo  volume  in  a  1  s  interval  prior  to  sprite  initiation.  We  note  that  (he  above  estimate  is  rough 
as  no  observations  are  yet  available  quantifying  the  meteor  flux  density  in  a  specific  region  above 
an  active  thunderstorm  producing  sprites.  Nevertheless*  this  estimate  suggests  a  higher  probabil¬ 
ity  (i.e.*  5%  using  the  methodology  described  above)  of  a  positive  correlation  between  sprite  and 
meteor  events  than  that  shown  in  previous  literature . 

As  a  summary*  a  comparison  of  experimental  and  numerical  results  presented  in  this  study 
demonstrates  that  the  descending  halo  structures  in  sprite-halo  events  recorded  by  high-speed  video 
cameras  are  the  optical  manifestation  of  pre-existing  plasma  irregularities  in  the  D-region  iono¬ 
sphere,  it  is  for  the  first  time  that  robust  evidence  of  plasma  irregularities  in  association  with 
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sprite  streamer  S initiation  is  presented  and  the  characteristic  size  of  these  plasma  irregularities  is 
quantified.  Thunderstorm  effects  and  meteor  events  are  suggested  to  be  two  possible  origins  of 
those  plasma  irregularities.  We  note  that  future  observations  focusing  on  halo  structures  and  mod’ 
eling  studies  are  necessary  to  understand  the  origin(s)  of  those  D-region  plasma  irregularities  in 
association  with  sprite  initiation. 

Methods 

Plasma  Fluid  Model.  In  our  plasma  fluid  model1*,  the  chemical  reactions  accounted  for  include 
the  electron  impact  ionization  of  N*  and  Gj,  the  electron  dissociative  attachment  to  O2 .  and  the 
electron  detachment  process  Q“+N*-re+N/),  Note  that  although  included,  the  effect  of  detach¬ 
ment  in  our  simulations  is  negligible  over  a  ^2  ms  time  scale.  Photoionization  processes  are 
included  using  the  three-group  SP;3  model 27 .  The  motion  of  chafed  species  is  simulated  by  solv¬ 
ing  the  drift-diffusion  equations  for  electrons  and  ions  coupled  with  the  Poisson's  equation The 
transport  equations  for  charged  species  are  solved  using  a  flux -corrected  transport  technique  that 
combines  ait  eighth-order  scheme  for  the  high-order  fluxes  and  a  donor  cell  scheme  for  the  low- 
order  fluxes33.  The  simulation  domain  extends  from  ground  up  to  95  km,  has  a  radius  of  95  tin, 
and  is  discretized  using  a  grid  with  a  spatial  resolution  of  ^120  m.  h  is  assumed  that  perfectly 
conducting  boundary  conditions  are  satisfied  on  all  boundaries  of  the  simulation  domain. 

Ambient  Electron  Density  Profile,  Tire  ambient  electron  density  profile  can  be  expressed  as:g: 

na{h)  =  1.43  x  [in"3]  f  D 


where  hf  [kmj  and  £  [km-3]  are  given  parameters  describing  reference  altitude  and  sharpness, 
respectively  We  assume  that  tf=85  km  and  £=0,5  km"1,  which  represent  a  typical  nighttime 
electron  density  profile.  The  positive  ion  density  iiiur  is  assumed  to  be  equal  to  the  election  density 
at  high  altitudes,  where  electron  density  is  higher  than  Iff*  nr3,  and  is  equal  to  JO8  m-*  at  lower 
altitudes.  Initial  ambient  negative  ion  density  is  then  calculated  based  on  charge  neutrality 

Lightning  Current  Waveform,  lire  waveform  for  lightning  current  moment  is  modeled  using  the 
following  formulation30: 

'«>=S(5)“4(i)'1  ® 

where  the  lightning  charge  moment  change  QIiq  is  chosen  to  be  300  C  km  and  hr 25  /rs  that 
corresponds  to  a  return  stroke  with  a  duration  of  ms. 

The  Image-Processing  Algorithm,  To  reconst  met  the  shape  and  size  of  the  plasma  irregularities, 
different  image  processing  algorithms  could  be  developed.  Here  we  provide  an  illustration  of  a 
simple  and  effective  algorithm  developed  as  part  of  the  present  work.  All  frames  from  the  begin¬ 
ning  of  the  halo  to  streamer  initiation  are  used  in  the  reconstruction  process.  First,  each  frame  is 
converted  to  a  8  bit  gray  scale  image.  Second,  the  intensity  in  each  image  is  linearly  normalized  to 
equalize  the  peak  intensity  in  different  images.  The  purpose  of  this  normalization  is  to  make  the 
brightness  of  the  halo  structures  in  different  images  that  corresponds  to  different  portions  of  the 
plasma  irregularities  comparable.  The  rationale  of  the  normalization  is  explained  as  follows.  The 
brightness  of  the  halo  structures  in  the  original  image  time  series  could  be  significantly  different 
front  each  other.  This  difference  is  mainly  due  to  the  temporal  variation  of  the  lightning  reduced 
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eleeiric  held  (see  e,g„  Figure  1  in  the  reference  [II]).  and  may  partially  stem  from  vertical  den¬ 
sity  variation  in  the  plasma  Irregularities,  In  other  words,  difference  in  the  brightness  of  the  halo 
structures  in  different  images  does  not  necessarily  indicate  that  the  density  in  the  body  of  a  plasma 
irregularity  varies  significantly  in  the  vertical  dimension.  As  shown  in  Figure  2d,  even  if  the  plasma 
irregularity  is  assumed  to  be  vertically  homogeneous*  temporal  variation  of  the  lightning  electric 
field  could  still  lead  to  significant  difference  in  the  brightness  of  the  halo  structure.  Since  our  aim 
is  to  reconstruct  the  shape  and  size,  rather  than  the  density  variation  of  ihe  plasma  irregularities, 
we  remove  the  effect  of  lightning  electric  field  variation  by  equalizing  the  peak  intensity  in  the 
sprite  images*  This  step  is  important  as  otherwise  some  pans  of  the  plasma  irregularities  cannot 
always  be  reconstructed  due  to  weak  brightness  when  compared  to  other  parts,  in  the  final  step*  we 
compare  the  normalized  images  pixel  by  pixel  retaining  only  the  maximum  intensity  in  each  pixel 
and  using  these  values  to  form  ihe  reconstructed  image.  Assuming  that  10  normalized  images  are 
used  in  the  reconstruction,  and  each  has  a  total  number  of  pixels  MxN,  For  each  pixel  denoted  by 
(mfir)*  there  are  10  different  intensity  values  from  those  10  normalized  images,  among 

which  ihe  maximum  one  »■)  is  chosen  to  be  the  intensity  in  the  pixel  (hi,  n)  of  the  recon¬ 

structed  image*  Using  this  algorithm*  the  halo  structures  captured  in  an  image  lime  series  can  be 
merged  in  a  single  image  to  reveal  the  shape  and  size  of  She  plasma  irregularities. 

L  Lay,  E,  H,  &  Shao,  X,  M.  High  temporal  and  spatial-resolution  detection  of  D-layer  fluctua¬ 
tions  by  using  time-domain  lightning  waveforms,  J.  Geophys ,  Rex,  116,  AO  1 317  (201 1). 

2,  Shao,  X,“M„  Lay,  E.  H.  &  Jacobson,  A.  R,  Reduction  of  electron  density  in  the  night-time 


lower  ionosphere  in  response  to  a  thunderstorm,  Nat.  Geosci.  6*  29-33  (2013). 

3.  Stenbaek-Nielsen,  H.  C,,  Kanmae,  T.,  McHarg*  M.  G,  &  Haaland,  R.  High-Speed  Observa¬ 
tions  of  Sprite  Streamers,  Sun*.  Geophys.  in  press  (2013). 

4.  Maihotra*  A,  &  Mathews*  j.  D,  Low-altitude  meteor  trail  echoes.  Geophys.  Rex ,  Lett.  36. 
L2!  106(2009), 

5.  Wilson*  C.  T,  R.  The  electric  field  of  a  thundercloud  and  some  of  its  effects.  Proc.  Phys ,  Sue. 
Loudon  37,  32D-37D  (1925). 

6.  Franz*  R,  C*  Nemzek*  R,  J,  &  W3 tickler,  J.  R,  Television  image  of  a  large  upward  electric 
discharge  above  a  thunderstorm  system.  Science  249*  48-5 1  (1990), 

7.  Pasko,  V,  P„  Inan*  U,  S,*  Roll,  T.  F.  &  Taranenko*  V*  N,  Sprites  produced  by  quasi-elect  rest  <Uic 
heating  and  ionization  in  the  lower  ionosphere,  J.  Geophys.  Res.  102*  4529-4561  (1997). 

8.  Luque,  A.  8l  Ebert,  U.  Emergence  of  sprite  streamers  from  screening-ionization  waves  in  the 
lower  ionosphere,  Nat.  Geosci.  2,  757-760  (2009), 

9.  Qin,  J.,  Celestin,  S*  &  Pasko*  V,  P,  On  the  inception  of  streamers  from  sprite  halo  events 
produced  by  lightning  dischatges  with  positive  and  negative  polarity.  J.  Geophys.  Res,  1 16, 
A063G5  (2011). 

10,  Cummer,  S.  A.  &  Lyons,  W*  A*  Implication  of  lightning  charge  moment  changes  for  sprite 
initiation.  7.  Geophys,  Rex.  1 10,  A04304  (2005), 


15 


16 


H.  Qin.  J.n  Ceieslin,  S.  &  Pasko,  V.  P.  Fo  miation  of  single  and  double- headed  streamers  in  sprite- 
halo  events.  Geophys.  Res.  Lett.  39,  L05SI0  (2012). 

12.  Lang,  T.  J.  et  al  Transient  luminous  events  above  two  mesoscale  convective  systems:  Charge 
moment  change  analysis.  J.  Geophys*  Res.  116,  A 10306  (2011), 

13.  Barrington-Leigh,  C.  P„  [nan,  Li.  S,  &  Stanley.  M,  Identification  of  sprites  and  elves  with 
intensified  video  and  broadband  array  photometry.  J.  Geophys.  Res.  106,  1741-1750(2001). 

14.  Pasko,  Vr  R,  Inan,  U.  S.  &  Bell,  T.  F.  Spatial  structure  of  sprites.  Geophys.  Res.  Lett ♦  25, 
2123-2126(1998). 

15.  Kanmae,  T„  Stenbaek-Nielsen,  H.  C,  McHarg,  M.  G,  &  Haaland,  R.  K,  Diameter-speed 
relation  of  sprite  streamers.  J.  Phys.  D-Appl.  Phys.  45,  275203  (2012). 

J6,  Liu,  N„  Kosar,  B„  Sadighi,  S..  Dwyer,  J.  R.  &  RassouJ,  H.  K.  Formation  of  Streamer  Dis¬ 
charges  from  an  Isolated  Ionization  Column  at  Subbreakdown  Conditions.  Phys.  Rev.  Lett. 
109,025002(2012). 

17.  Kosar,  B.  C„  Liu,  N.  &  Rassoul,  li.  K.  Luminosity  and  propagation  characteristics  of  sprite 
streamers  initiated  from  small  ionospheric  disturbances  at  sub- breakdown  conditions.  /  Geo¬ 
phys.  Res.  117,  A08328  (2012). 

18.  Qin,  J„  Celesta n,  S,  &  Pasko,  Vr  P.  Dependence  of  positive  and  negative  sprite  morphology 
on  lightning  characteristics  and  upper  atmospheric  ambient  conditions.  J.  Geophys.  Res.  1 18, 
2623-2638(2013), 


17 


19.  Marshall,  R.  A.  Ait  improved  model  oflhe  lightning  electromagnetic  field  interaction  with  the 
D-region  ionosphere.  J.  Geophys.  Res.  I17r  A 03 3 1 6  (20 ] 2). 

20.  Cordillo- Vazquez.  F,  J.  Air  plasma  kinetics  under  the  influence  of  sprites.  J.  Phys ,  D:  AppL 
Phys.  41,234016(2008). 

2L  Stenbaek-Nielsen,  H.  C,  Moudry,  D.  R.,  Wescott,  E.  M.,  Sent  mao,  D,  D,  &  Sabbas,  F.  T,  S, 
Sprites  and  possible  mesospheric  effects.  Geophys.  Res.  Lett.  27,  3829-3832(2000). 

22,  Sentman,  D.  D„  Stenbaek-Nielsen,  H.  C„  McHarg,  M.  G.  &  Morrill,  J.  S,  Plasma  chemistry 
of  sprite  streamers.  J.  Geophys.  Res.  113,  Dll  112  (2008). 

23,  Haldoupis,  C.,  Cohen,  M„  Cotts,  B.,  Arnone,  E.  &  Inan,  U.  Long- lasting  D-region  ionospheric 
modi  heat  ions,  caused  by  intense  lightning  in  association  with  elve  and  sprite  pairs.  Geophys. 
Res.  Lett.  39,  LI6801  (2012). 

24,  Suszcynsky,  D.  M  et  al.  Video  and  photometric  observations  of  a  sprite  in  coincidence  with  a 
meteor-triggered  jet  event,  J.  Geophys ♦  Res.  104.  3 136 1 -3 1 368  (1999), 

25,  Verbeeck,  C,  Meteoroid  flux  densities  from  radio  observations:  the  latest  developments 
(2009).  International  Meteor  Conference,  The  International  Meteor  Organization,  Porec,  Croa¬ 
tia.  24-27  Sept, 

26,  Wescott,  E.  M,  et  al.  Triangulation  of  sprites,  associated  halos  and  their  possible  relation  to 
causative  lightning  and  micrometeors.  J.  Geophys.  Res.  106,  10467-10478  (2001), 


27.  Bourdon.  A*  et  ctl.  Efficient  models  for  pbotoionization  produced  by  non-thermal  gas  dis- 
charges  in  air  based  on  radiative  transfer  and  the  Helmholtz  equations.  Plasma  Sources  ScL 
Technol.  16,  656-678  (2007), 

28.  Zalesak*  S,  T.  Fully  multidimensional  flux -corrected  transport  algorithms  lor  fluids.  X  Com- 
pitL  Rhys,  31,  335-362  (1979). 

29,  Wail,  J,  R.  &  Spies,  K.  P.  Characteristics  of  the  Earth-ionosphere  waveguide  for  VLF  radio 
ivmw,  Tech  note  300  (NaiionaJ  Bureau  of  Standards,  Boulder,  Colorado,  1964), 

30,  Clio,  M,  &  Rycroft,  M.  J.  Computer  simulation  of  the  electric  field  structure  and  optical 
emission  from  cloud- 1  op  to  the  ionosphere.  X  Almas,  Solar  Ten:  Phys,  60+  871-888  (1998), 

Acknowledgements  This  research  was  supported  by  the  Defense  Advanced  Research  Projects  Agency 
and  by  the  Physical  and  Dynamic  Meteorology  and  the  Agronomy  Programs  of  National  Science  Founda¬ 
tion.  We  (hank  I.  Mathews  for  useful  discussion  about  possible  relation  between  sprites  atid  meteors. 

Author  contributions  J.Q.  drafted  the  manuscript,  developed  the  plasma  fluid  model*  interpreted  simu¬ 
lation  results  and  compared  them  with  the  observations.  V.RP  supervised  the  project.  M.G.M  and  H.C.S. 
operated  (he  sprite  observations,  derived  the  altitude  range  of  the  images.  All  authors  contributed  to  tlte 
discussion  of  lire  results  and  the  preparation  of  the  manuscript. 

Competing  financial  interests  The  authors  declare  no  competing  financial  interests. 

Correspondence  Correspondence  and  requests  for  materials  should  he  addressed  to  Jianqi  Qin  (jian- 
qiqin@psu.edu). 


19 


-20 


Radial  position  (km) 


-20 


■ .  s  '*  - 


|  fd)  image  series 


Sprite  streams* 


*  J**l 

J*4| 


0.0  05  1.0  ),S  1.375 

Time  (ms) 


Figure  It  A  sprite- halo  event  observed  on  August  27,  2009.  The  event  starts  at 
09:15:23  085  570  UT  (t=U  ms).  The  frame  rate  of  the  camera  is  10,000  fps,  Images  (a,  l>) 
respectively  show  the  halo  emissions  at  £=0,750  ms  and  1,125  ms,  between  which  the  halo 
structure  descends  with  an  average  downward  speed  of  6,4  x  Ilf  m/s.  Image  (c)  shows  the 
sprite  streamer  initiated  from  the  descending  spatial  structure.  Panel  (d)  shows  a  1,875-tns 
image  time  series  of  the  halo  structure  and  streamer  development. 
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Figure  2:  Descending  halo  structure  observed  in  a  halo  modeling,  (a)  Ambient 
electron  density  and  the  plasma  irregularity  In  the  simulation  domain .  (b,  c)  Halo  etnissiojis 
from  the  first  positive  band  system  of  Na  at  f— 0.7  ins  and  1.2  ms,  respectively,  between 
which  Hie  bright  portion  descends  downward  at  average  speed  of  7.0x10°  m/s.  (d)  A  1,4-tns 
image  time  series  of  the  modeled  halo  structure. 
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Reconstructed  ptesma  irregularity 
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Figure  3:  Reconstructed  plasma  irregularities,  (a)  Tire  plasma  irregularity  recon¬ 
structed  using  image  series  from  t=Q  ms  to  1.4  ms  corrcspoiuling  to  the  model  event  sliown 
in  Figure  2.  (b)  The  D-rcgion  plasma  irregularity  reconstructed  tising  image  series  of  the 
sprite  halo  event  shown  in  Figure  !♦ 
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Figure  4:  Halo  events  with  and  without  plasma  irregularities,  (a,  b)  A  frame  of 
the  sprite  halo  event  observed  on  20  July  2012  and  the  reconstructed  Image  showing  five 
plasma  irregularities*  (cT  d)  A  frame  of  the  halo  event  observed  on  14  July  2010  and  the 
reconstructed  image  that  shows  only  large-scale  halo  emissions. 
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